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Attenuation of the singlet quark-anti-quark bound state passing through a random external SU(7V) 
color field at high energy is studied within the eikonal framework as a model for the absorption 
of a quarkonium by nuclear target. For a thick target of length L the penetration probability 
is suppressed non-exponentially and behaves asymptotically cx with indicating the 

equi-partitioning in the color space. An approximate rational expression for the penetration of the 
small Gaussian bound state is found which allows the explicit L and A*' dependences, and shows the 
attenuation is in the powers of L although the approach to the asymptote is rather slow. 



PACS numbers: ll.80.La, 13.85.-t, 24.85. +p 
I. INTRODUCTION 

Color dipole is a simplest physical (color singlet) state 
of QCD, and hence may be used as a tool to study more 
complex objects in QCD. In Nature it is realized as heavy 
quarkonia (cc or bb) , a virtual qq pair in deeply inelastic 
scattering seen in the target rest frame, and so on. In the 
context of the nucleus-nucleus (AB) collision physics, the 
modification of the quarkonium productions in the QCD 
media has been extensively studied both experimentally 
and theoretically ^ since the proposal as a signal of the 
QCD plasma formation H. 

In the previous paper [g| we studied the penetration 
of a high energy QQ bound state traversing a nuclear 
target, which is modeled as soft random classical fluc- 
tuations of SU(2) gauge fields. Working in the eikonal 
framework we showed that the penetration probability 
of the color singlet bound state attenuates asymptoti- 
cally in 1/L with L the target thickness |^,|^ , not in the 
exponential form. This is essentially the same effect as 
the so-called superpenetration, conjectured theoretically 
for the ultra-relativistic positronium passing through a 
metalhc foil |,|,|. 

The aim of this paper is to extend our formulation 
to general SU(Af) case. We explicitly show the 1/N'^L 
suppression of the penetration probability of the singlet 
bound state in the asymptotic region. An analytic, ap- 
proximate formula which well reproduces the N- and L- 
dependences is given for the Gaussian bound state and 
shows that the penetration is power-suppressed by 1/L 
although the approach to the asymptotic form is rather 
slow. 

Within the eikonal approximation we deal with the 
multiple scatterings of the quark and anti-quark with 
the phase factors accumulated along the classical straight 
paths with fixed transverse positions. The bound state 
wavefunction described in terms of these phases should 
generally contain various excited states of the pair from 
the initial state. In this sense our treatment of the 
pair propagation on the background field includes the 
quantum-mechanical coherence. The random phase dis- 
arranges the wavefunction to reduce the overlap with the 
initial state on average. 



In Sec. H the basic eikonal formulas for SU(iV) the- 
ories are given, and Sec. Ill defines the target average 
following Ref. . Our main result on the survival prob- 
ability of the singlet bound state is presented in Sec. IV. 
Sees. ^ and [V^ are devoted to discussions and summary. 
In Appendix |A| the target averaging of the pair eikonal 
factors is briefly explained. Our framework allows an an- 
alytic solution to the SU(2) model with the finite gauge 
fluctuation amplitude of Ref. ||^, which is given in Ap- 
pendix pi 



II. EIKONAL APPROXIMATION FOR PAIR 
PROPAGATION 

Let us consider the quark-anti-quark bound state pro- 
jected onto a large target nucleus of thickness L at very 
high energies. We choose the target rest frame in which 
the projectile bound state travels with ultra-relativistic 
momentum P. In such a high energy process the de- 
flections of the constituents of the bound state can be 
ignored and their transverse positions are frozen during 
the interaction with the target. The eikonal approxima- 
tion for the propagation of these constituents is our basic 
framework to describe the interactions of the bound state 
with the soft classical gluon fields created by fluctuations 
in the nucleons in the target nucleus |10Hlq|. 

The bound state is described by a wavefunction 
tf{r, z)^ where r = x — y is the relative transverse po- 
sition of the quark at x and anti-quark y, z is the frac- 
tion of the longitudinal momentum carried by the quark, 
and ^ is the color wavefunction of the pair. Except for 
the vicinity of the endpoints, z = and 1, the longitu- 
dinal momenta of the constituents are large enough that 
the quark and anti-quark are regarded as moving along 
the classical, straight-line trajectories. The interaction 
with the external gauge flelds is described by the usual 
eikonal quark-gluon vertex igv^A^(x) with = p^/Ep 
and = (t° are the SU(A^) generators normal- 

ized as Trt°t^ — ^5°-^). Traveling through the target 
nucleus, the quark acquires the eikonal phase, or Wilson 
line, VF(x) accumulated along the path as 
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\g J"^ dx+A" (x,an+)t" 



(1) 



where we have taken the hmit — (1,0,0,1) and in- 
troduced the notations a;* — -^{t ± z) and = 

--^(v4° zL A^). This is a matrix in the color space and 

V indicates the path ordering of the product. The cor- 
responding eikonal factor for the anti-quark is simply 

W^y;A) = 7'e'''^+~ since it belongs to 

the color representation conjugate to the quarks. After 
the interaction with the target nucleus, hence, the quark- 
anti-quark pair gets the total eikonal phase 



U{^,y;A)^Wi^;A)wHy;A). 



(2) 



The scattering amplitude from the initial state 
(/>i(x, y, z) = ipi{r,z)^i to the final state 0f(x, y, z) = 
e~''*^^(y5t(r, z)^f with the CM transverse momentum Q 
and R = ^(x -I- y), is expressed as 



($f0f I f/| 



DAdyLdjdz^ [A]e'-^^ip*^ (r, z)i\ 
x;7(x,y;A)$i[A]^i(r,z)6, (3) 



where ^n^) is the initial (final) state of the target which 
we treat as a functional of soft external fields, A. 

In studying the penetration probability of the initial 
bound state, only the out-going state of the pair is of 
relevance. Taking the closure of the final target states, 
we have 



(<^i|m|0f)((/.f|c/|(/.i) = (<i>i| (</)i|f/t|(/.f)(^f|f/|(/.i) |<i>i) 

DA {4>,\U^[Am{^,\U[Am. (4) 



In this form we first calculate the scattering probability of 
the pair on a certain background field configuration, and 
then take average over the gauge configuration distribu- 
tion It is more convenient, however, to average 
over the A configurations before taking the matrix ele- 
ment with respect to the pair bound state, because the 
A fields appear only in the eikonal factors which are in- 
dependent of the in- and out-states. More explicitly it is 
written as 



(0i|C/t|0f)(0f|{7|0i) ^ / dxdydzdMyd. 



„iQ(R-R.) 



x.^f(r, z)(^f(r,z)Pfi(x,y;x,y)(pi(r,z)<^i(r,z) , (5) 

where x, y and z are the coordinates in the complex 
conjugate amplitude, and we introduced a shorthand no- 
tation for the product of the eikonal factors. 



Pfi(x,y;x,y) = 4^f/(x,y)ei (^J C/(x, y)^) 



(6) 



Pfi includes all the dependence of the probability on the 
background gauge field, and the probability is obtained 
as an integral with the kernel Pfi. 




nuclear target 

FIG. 1. Schematic picture of the high energy QQ pair pass- 
ing through a nuclear target. 



A. eikonal factors in color irreducible basis 

The expressions for the eikonal factors are much simpli- 
fied in the color basis of the irreducible representations. 
Since the color singlet and adjoint states are constructed 



as 
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7^ L.^=l and ^ = V2E„=i N)i?,01, re- 

spectively, out of the quark |i) and the anti-quark 
the elements of the eikonal factor with the definite color 
yield 

C/,,(x,y) - j^TT{W{^)WHy)), 



U:j^,y) = ^ j^TviWi^)WHyr), 

U,liK,y) = y^Tr(VF(x)rM/t(y)), 
UtA^,y) = 2Tiit'W{^rwHy)). 



(7) 



These eikonal factors show the crucial features of the 
interactions of the pair. Color Transparency p^ , p^ : An 
elementary color-singlet dipole with a vanishing size, x = 
y, does not interact with the gauge fields at all. Actually 
the quark and the anti-quark acquire the same eikonal 
phase, and hence C/gs = 1 and Uas — 0, irrespective of 
the gauge fields. Color Rotation: An elementary color- 
adjoint dipole changes its color orientation as well as the 
CM motion through the interaction with the gauge fields. 

This is an x matrix in the color space. In Ref. |^ 
we explicitly dealt with this full matrix for SU(2) case, 
and confirmed that the color rotation within the adjoint 
representation is secondary for the dynamics of the pene- 
tration of the color singlet state. Furthermore the orien- 
tation of the color adjoint vector has almost no meaning 
in the color gauge theory. Hence we will concentrate 
on the transition processes between the irreducible color 
representations of the singlet and the adjoint. 

Summing up the color orientations, a = 1, • • • , iV^ — 1, 
we have the eikonal factors squared for the singlet-singlet 
(Pss), singlet -adjoint (Pas and Pga), adjoint -adjoint (Paa) 
processes, which are compactly expressed with the traces 
of the products of the Wilson line factors as follows: 

Pss = [/ss(x,y)C/,;(x,y) ^ ^Oi(x,y,x,y), (8) 
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Pas^Ef^as(x,y)(7r(x,y) 



^02(x,y,x,y)-^0i(x,y,x,y), (9) 



F,a^^C/,1(x,y)C/r(x,y) 



l02(x,x,y,y)-^0i(x,y,x,y), (10) 



P,,^^C/inx,y)C/,^r(x,y) 



= C'i(x,x,y,y) - — 02(x,y,x,y) 

-l02(x,x,y,y) + ^0i(x,y,x,y), (11) 

where we defined 

Oi(x,y,x,y) = Tr[PF(x)PFt(y)] TY[W{y)WH5c)], 
02(x,y,x,y) EE Tr[W^(x)W^t(y)^(y)^t(5^)]. (12) 



B. cross sections 

For the sake of simplicity, we assume here the nu- 
clear target with infinite transverse extent hypothetically. 
The survival probability for a color singlet initial state, 
(/3o(r, z)£,s, is expressed as 

S = dxdydzdxdydz ^^(R - R) 



x.^o(r, z)^o(r, 2)Pss(x, y; x, y) 



drdrpo(r)po(r)P5s(r,r,0) , 



(13) 



where A is the normalization transverse area and the 
S function appears after the integration over Q. Here 
we also introduced /cio(r) = JdzipQ{r,z). Note that 
Pg(x, y;x,y) = Pfi(r,f,R — R) due to the assumption of 
the translation invariance in the transverse plane of the 
target. Then the normalization area A is canceled out 
correctly. 

At this point it may be useful to give the expressions 
for the total and diffraction cross sections as well as the 
elastic cross section of the singlet state. The total cross 
section is obtained by taking the closure of all the possible 
pair final states as 



cr^* = J diidydz(pl{r,z) 



x-Tr[(C/(x,y)-l)(C/t(x,y)-l)] 



A / drpo(r) (2 - 2ReC/s.(r)) 



(14) 



where Uss (x, y) = Uss (r) , and the sufRx A indicates that 
this is the contribution from the transverse area A. It 
is clear that the factor 2=1+1 represents the incoherent 
sum of the contributions of the incident and scattered 
waves, and the second term expresses the interference. 

The diffraction cross section, which is defined here as 
the processes where the pair remains color-singlet with- 
out the net color exchange with the target, but can be ex- 
cited by the interaction, is obtained by adding the color- 
singlet condition on the out-going pair: 



.diff 



dxdydz(pl{r, z) 



x(;7..(x,y)-l)([/*(x,y)-l) 
A [ drpo(r) (i^(r, r, 0) - 2ReC4l(r) + l) . (15) 



Finally the differential elastic cross section is 



dai 



dQ/(27r)2 

= J dxdydzdxdydzipl{r,z)ipl{r,z)e'^'^^~'^'^ 

x(C/,,(x,y)-l)(C/ls(i,y)-l) 
= aJ d(R-R)drdfpo(r)/5o(f)e*'5(^-^) 

X f7^(r,f,R-R)-[7:r(r)-C4(f) + l) . (16) 



We note here that the inelastic cross section defined by 
CT^ = (7^* — cr^, satisfies the usual relation 

<j'J=aJ drrfrpo(r)po(r) (l - 7^(r, f , 0)) . 

= A{1-S). (17) 



III. MODEL OF THE TARGET 

We are considering the target as a source of the back- 
ground gauge fields created as the fluctuations in the nu- 
cleons in the thick target of length L. We divide the path 
of the pair into n different small zones. The size I (—L/n) 
of the zone may be taken as the correlation length of the 
fluctuations, which should be of the order of, or smaller 
than the nucleon size in the target rest frame. We assume 
here that the flelds in the different zones are uncorrelated 
and then we take the target average in each zone inde- 
pendently ||l^,|l3|,||. The eikonal factor is now factorized 
into a product of the contribution from each zone: 



W{x.) = M^(x)(„) = VF„(x) • • • M^2(x)W^i(x). 



(18) 



We indicate here, and hereafter if necessary, by the lower 
suffix (n) the fact that the quantity is a product of n 
factors from the independent zones in the target. In 
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the transverse direction, we define the distribution of the 
color field fiuctuations by 



A-{^)A\y) = 5'''C{^-Y) , 



(19) 



where A" = gAx+A"- with Ax+ = ^21. Note that the 
transverse correlation length should be the same order of 
I for consistency. 

We follow Ref. about the field averaging and the 
details are presented in Appendix We expect that the 
fiuctuating field A in each zone is weak and expand the 
eikonal amplitude up to 0{g^). For example MS, the 



eikonal factor appearing in the total cross section (14) 
averages 



C/ss(x,y) 



1 



Tr(W^(x)W^t(y)) 



N' 

[1 - CMrT 



(») 

-nCFv{r) 



(20) 



where Cf — 



2N 



and v{y) = C(0) - C(r) > the 
subtracted correlation function. 

In the diffraction cross section, which yields the con- 
dition X = X and y = y for the transverse positions in 
the amplitudes and the conjugate one, the eikonal factors 
Pg(x, y; X, y) are reduced to 

P..(x,y;x,y) = \Uss{^,y)\\ 

Pas(x,y;x,y) = Psa(x,y;x,y) = 1 - |C/ss(x,y)p, 

Paa(x,y;x,y) =iV2_2+|[/,,(x,y)|2. (21) 



The target average results in a compact form 
P..(x, y; X, y) = — [1 + {N^ - 1)(1 - Nv{v)r] 



1 
1 

7V2 



1 



(22) 



and those for other channels follows from the above rela- 
tion (^. Eq. ( p2| ) is regarded as the probability that the 
singlet pair of size r remains singlet. It shows that the 
equal-probability (1/7V^) distribution in the color space 
is achieved exponentially in L = nl for the quark pair 
with a definite size r B . The result of 



2Re[/,,(r) + 1 
+ 1 - 2e~"'^^"("") + 



Pss(r,r,0 
1 

iV2 



-nNv(r) 



(23) 



is already given in Ref. g] . 

With this model of the target we obtain the recursion 
formula for the average of the eikonal factor Pfi(x, y; x, y) 
(see Appendix 0), which is easily solved as 



P P 

aa -'as 
P P 



(n) 



P P 

aa ^ af 

P P 



(n-l) 



~1 
1 



(24) 



The fact that the equation is closed in Pa (f , i = s, a) only, 
shows that the transition between the singlet and adjoint 
states of the pair is a stochastic process during the travel 
through the target. That is, the non- vanishing contri- 
butions in the color averaging come from those where 
the color states of the pair in the amplitude and in the 
conjugate amplitude are the same in each step. 
The explicit matrix elements of U are 

Z^aa = l-CF[t;(x-x) + t;(y-y)] 

+ - S) + v{y - y) + v{x - y) 

+v{5c - y) - 2v(x - y) - 2w(y - x)], 
Ui,s = Cf[w(x - y) + v{y - x) - w(x - x) - v{y - y)] , 
1 



_ I 1 

Uss^l- Cf[v{x -y) + u(S - y)]. 



(25) 



As we summed over the orientations of the adjoint vec- 
tors, Uas is — 1 times as large as Usa 

Since U is a 2 x 2 matrix, the calculation of U"^ is a sim- 
ple task once the explicit form is given. For a thick target 
(large n), U"' may be sharply peaked at x = y = x = y, 
where U is unity, and therefore we can approximate it 
with a Gaussian in the sense of the steepest descent: The 
explicit form of U we need in this approximation is 



Ua 



1 



+ (iV^-3)(r-r)2 + (r + r)2], 
Uas = iCpKr ■ f , 

7/ ^ - 

Z^sa = ^r.r, 



(26) 



with i;(r) ~ Kr^ and k oc {V j_A"-{0)y . We note again 
that the color transparency effects are properly included 
in the fact that Uss = 1 and Uns = ^sa = for r and/or 
f = 0. From this expression for U we find explicitly that 
Pfi is a function of r, f , R — R as noticed in the previous 
section based on the translation invariance of the target. 
The R — R dependence appearing in U^a should lead 
to the diffusion of the CM transverse momentum of the 
surviving bound state. 



IV. SINGLET BOUND STATE PENETRATION 

In the thin target limit, the interaction should be de- 
scribed with 71 = 1, that is, single zone on the path. 
Noting that in the thin nuclear target the inelastic cross 
section is a sum of those of the individual nucleons in 
the target cr^ = ^Zriocr'" with no nuclear number den- 
sity, we can directly relate the expectation value of the 
correlation function ^(r) to the inelastic cross section of 
the singlet bound state with a nucleon through 



4 



(Alno) d'" = J dRdrdfpn{r)po{f ){l-Uss) 

= ^ j drdfpQ{r)po{r)CF{v{r) + v{r)) 
= 2ACf{v) (27) 

with (v) = J (irpo(r)ti(r). 

Converting (v) with the physical quantity ct'", we find 
that the survival probability ( p^ for a small L ^ nl has 
the usual form S = 1 - a"^noL = 1 - L/Lin — e^^^^^'^ 
with 1/Lin = cr'"no = 2Cf{v)1~^. 

For a thick target the eikonal factor to calculate the 
survival probability of the singlet bound state is ex- 
pressed as 



f'ss(r,r,0) = /o;ie 



-nln(l/Ai) 



/0;2e 



-riln(l/A2) 



(28) 



where /o;i,2 are the singlet components squared of the 
eigenvectors associated with the eigenvalues Ai_2, respec- 
tively, of the matrix U. For a thick nuclear target to 
which the approximation ( p6[ ) applies, the explicit forms 
for these factors are found to be 



1 



2r2 + (Ar2 _ 2)? 



/o;i,2 2 ^ 2N[Art + iVV^ - r2_]i/2 



(29) 



Ai,: 



K 

8N 



2(7V^ - 2)rt +N^rl 



4^2 ^2 11/2 



(30) 



with r| = (r ± r)^. The eigenvalues satisfy A2 < Ai < 1, 
and Ai = 1 when r = r, which corresponds to the unit 
eigenvalue of the probability matrix ( |A^ ) in the diffrac- 
tion. With N = 2 this reproduces the result of Ref. j|] . 

In the asymptotic region the larger eigenvalue Ai con- 
trols the behavior of the probability. Noting that in the 
exponent 

ln(l/Ai) ~ KCi.r2/2 , 



the asymptotic behavior of the survival probability is ob- 
tained as 



drdrpo(r)/Oo(r)/o;ie" 



N 
1 



A 



7V2 nnCp ' 



(32) 



where A = 27r j drpQ{r). In SU(2) the fundamental rep- 
resentation is pseudo-real and Eqs. (^9|) and (^) become 



symmetric in with iV = 2, which multiplies an 

extra factor of 2 on RHS ||^ . The most important fact we 
observe here is that the integration over r with the weight 
po(r) changes the superficial exponential dependence on 
n (X L into the power dependence, l/n oc 1/L. 

For a small bound state where {v) ~ K(r^) is valid even 
for a thin target, the single parameter n determines the 




FIG. 2. Survival probability of the color singlet Gaussian 
bound state passin g th rough SU(3) random background fields. 
The analytic form ( p5| ) [solid] and the simple exponential form 
[dashed] are compared with the numerical result [cross]. The 
first term of Eq. (bsf) is also shown in dotted line. 



behavior of the survival probability from the thin to the 
thick target region. With ii„ the asymptotic form (|32" 
is rewritten as 



5: 



A 1 

7V2 2^ 



-Lr^_/4Li„{r^) _ 



= 2 



A(r2) U 



7V2 



(33) 



In this form the N dependence appears in the denomina- 
tor only, which comes from the equi-partitioning in the 
color space. 

It would be suggestive to show an explicit calculation 
with the Gaussian wavefunction po{r) — ^^^^^^ Z^"" ^ 

and using the form w(r) = Kr^. For the Gaussian wave 
function A(r^) = 1. The survival probability scales as 
S = S (L/Lin) since the exponent in the eikonal kernel is 
quadratic in r and r and we can rescale the integration 
variable to r/(r2)i/2 and f/(r2)i/2 (see Eqs. (H), (|o|)), 
(31) which results in a prefactor, nK(r^), in the exponent of 



the eikonal kernel. 

In U(l) case this simply yields 



S 



1 



l + 2nK(r2) 1 + L/Lu 



(34) 



where Lin = 1/{2kI^'^{P)). Eq. (||) smoothly interpo- 
lates between the usual absorption rate 1 — L/L\-a for 
small L and Li^/L for large L. 

In SU(7V) we know the asymptotic form both in the 
large and small L regions from the above results in this 
section. Furthermore, the integral (|l^) with ( p8| ) in the 
region except for r ~ f should behave as for large L 
since the exponent of the eikonal kernel is proportional 
to the target width L. From this consideration and mo- 
tivated by the U(l) result we speculate an analytic form 
for S as follows: 



S{L/Lin) 



1 



-2 



1 



(35) 



5 



This reproduces the correct small and large L behavior 
of 5. 

We compare this form (p5) with the result of the nu- 
merical integration in Fig. ^. We find that the sim- 
ple analytic form ( |35| ) approximates the numerical re- 
sult very well. We confirmed this for = 3 ^ 10 and 
L/Liri = ~ 10. The fit in SU(2) is not as precise as in 
SU(iV > 3), due to the extra factor of 2 in the asymp- 
totic region for SU(2), but remains very much reasonable 
in the region of L/Li^ = Q ^ 2. 

In summary the survival probability of the singlet 
bound state decays in 1/L in the asymptotically large 
L. However, the approach to this asymptotic behavior 
is rather slow and the second term in Eq. ( p5| ) equally 
contributes in the intermediate region of L/Lin- For the 
larger N, it can even dominate over the first term. Nev- 
ertheless, the suppression of the survival probability is 
given in the power of 1/L, as typically shown in Eq. (35), 
not in the exponential one. 



V. DISCUSSION 

Within the same approximation of (v) = K(r^) and 
with the use of the Gaussian wavefunction, the asymp- 
totic forms for the total and diffraction cross sections are 
easily found as 



2-2e" 




(36) 



and 



difl 



A 



^'^^o(r) ( + l 



A 



2Cf i^-^^ 

N 

1 

iV2 + 



2Cf 



2e 
1 



TV 1 + J^-^ 

2Cf Lin 



1 + i-^ 

^ 2 Lin 



(37) 



We note that the above results show that cr^* — Aj^ 
and cr^'^ = 0{{^Y) for smaU L. The factor ^-^^ m 
the exponent of the integrand, or more generally 



w(r) w(r^) 



a no = a (r )no, 



(38) 



might be interpreted as a size-dependent cross section of 
the rigid dipole of size r. This fact that the smaller com- 
ponent is less interactive in passing through the target 
is usually referred to as color transparency. The integra- 
tion over the incoherent contributions of the color dipoles 



with size r leads to the power behavior in L for the to- 
tal and diffractive cross sections Q. In the total cross 
section (^6|) the interference term between the initial and 
scattered waves decays in 1/L, approaching the black- 
disk limit (73°^ — 2A. The diffraction cross section also 
approaches the incoherent sum of the incident and scat- 
tered wave contributions, A{1 + the latter of which 
is suppressed by the equi-partitioning factor in the 
color space. 

The asymptotic L^^-dependence of the survival prob- 
ability looks apparently the same as the one seen in the 
total cross section. The exponent of the integrand in 
Eq. (28), however, no longer has the simple meaning of 
the cross section tT'"(r^) of the dipole of size r, but de- 
pends on both of r and r in the amplitudes. Actually 
there are two exponential factors in the eikonal kernel; 
Ai which controls the asymptotic behavior of S, weights 
the diagonal components (r ~ f) of the amplitudes, and 
the other expresses the color diffusion in the color space. 

The question about the energy scale where the change 
of the survival probability from the usual exponential 
form to the power form occurs should be addressed. For 
this purpose one would have to includes the longitudinal 
momentum transfers and the binding energies of the res- 
onances. We leave this extension of the framework for 
future study. 

Another energy dependence will come from the split- 
ting of the quark to a quark and a gluon or evolution of 
the the quark pair wavefunction [EqjSl . Such an evolution 
may be also necessary to our model at very high energies. 

In the infinite momentum frame of the target, on the 
other hand, this branching should be included in the tar- 
get wavefunction. In McLcrran-Venugopalan model [T^ ] 
the valence charge distribution is treated as the random 
variable in the longitudinal as well as the transverse di- 
rections of the target. The generated classical fields has 
the subtracted correlation v{r) = Kr^ ln(l/r^A^) for suf- 
ficiently small ^ with A the infrared cutoff, re- 
flecting the massless nature of the fields. The inclusion 
of the quantum corrections to this model has been inten- 
sively studied at very high energies (see e.g., |l|]). This 
energy dependence might be economically mimicked by 
using the energy-dependent correlation function v(r) in 
our framework. 

As for the J/iIj production in heavy ion collisions (HIC) 
we need the initial quark pair wavefunction whose size 
is typically of the order l/2mc. At very high energy 
there is no time for the quark pair to expand, and the 
quarks stay close to each other during two nuclei going 
through. In this situation, the color fields in the nuclei 
will affect the penetration of the singlet state only weakly 
due to the color transparency. From the same reason, the 
transitions between the singlet and the octet states will 
be suppressed as seen in Eq. ( p^ ) ]l9| . 

One can argue based on the eikonal factor (|6|) that 
the CM momentum of the octet pair will diffuse stronger 
than that of the singlet one. The difference in the CM 
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momentum dependence between the singlet and adjoint 
states may give a hint on the production mechanism of 
the 3/tjj in the HIC experiments. 

Finally we notice that the J/ip production rate should 
given by the square of the total amplitude. Here we con- 
centrated only on the interaction of the QQ pair with the 
classical color fields in the target nucleus, ignoring other 
effects; for instance, the action of the background gauge 
fields on the quark pair may induce the gluon emissions, 
which will be important for the energy-loss and the color 
neutralization of the pair. We also leave these effects for 
future study. 



VI. SUMMARY 

In this paper we studied the absorption of the high- 
energy quark-anti-quark bound state in the random, 
classical, background SU(Af) gauge fields. Using the 
eikonal approximation for the quark and anti-quark prop- 
agations, we found that the probability is suppressed by 
the inverse power of the thickness L of the target for gen- 
eral SU(iV). The asymptotic form, was given, but 
its approach is rather slow in L. For the small Gaussian 
bound state, we interpolated between the known asymp- 
totic forms in the small and large L regions by the ra- 
tional form (^5|), and found that it gives a remarkably 
precise description for the numerical result in the whole 
i/iin region. At very high energies the survival proba- 
bility cannot have the exponential suppression behavior. 

In the inclusive cases such as the total and the diffrac- 
tive cross sections, the cross section of the bound state 
may be regarded as a sum of the dipole components of 
definite size r. However in more exclusive cases, e.g., the 
penetration probability S*, this simple interpretation does 
not apply. 

The survival probability of the J/i/j has an importance 
in the HIC physics as a QGP sensitive probe. In order to 
investigate the issue of the i/ijj yield in HIC more real- 
istically, we must certainly extend our model to include 
the processes of the heavy-pair production, the interac- 
tions with other hard constituents in the hot matter, and 
so on. 
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APPENDIX A: TARGET AVERAGES 

We follow Ref. in the target average. To perform 
the target average for the last zone we expand the eikonal 
factor or the Wilson line up to 0{g^) 



M/(x)(„) = VK„(x)iy(x)(„_i; 
1 + zA"(x, n)e - ^C{Q) ] I^(x)(„_i). (Al) 



The target average is defined by 



^"^(x, n)A''(y, n) = (5°''C(x - y) 



(A2) 



and we utilize an elementary relation = \(5ii5jk — 

jj^ijSki) for the fundamental representation. In the 
above expression we have averaged the 0{g^) term as we 
are working only up to this order. With this definition 
for the target average, we obtain 



t/s.(x,y) = -Tr(I^(x)W-t(y))^^^ 



-Tr(T4^(x)M^t(y))^^_^^[l _ CMr)] 



[1 -CFw(r)]" ~e-"'='^"W. 



(A3) 



The target average for the eikonal kernel for the sin- 
glet survival probability Pss is expressed as a trace of 
four Wilson line factors. First let us perform the target 
average of Oi and O2 given in Eq. (12). By taking the 
average for the last zone we find 



O2 



(x,y,x,y) 



(n) 



(x,y,x,y), (A4) 



(n-l) 



where 



a = 1 - Cf(u(x - y) + t;(x - y)) 

/3 = ^(x - y) + i;(x - y) - u(x - x) - i;(y - y) 

a' = 1 - Cf{v{-k - x) + i;(y - y)) 

/3' = t;(x - y) + w(x - y) - t;(x - y) - ^(x - y). (A5) 

Making the physical combinations, Pfi(x, y; x, y), 
given in Eq.(0), we find a closed relation 



(") 



2N ' 
— d 



(A6) 



(n-l) 



which is given in Eqs. (|24|), (pSj). The consistency may 
be checked by taking the average with rt = 1 on LHS and 
noting that Paa(o) = iV^ - 1, Pas(o) = Psa(o) = 0, and 
P,s(o) = 1 on RHS. 
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The target average of the eikonal factors Pa (x, y; x, y) , 
which appears in the diffraction case with f = i = s, are 
easier to do. Actually Pss(x, y: x, y) is a single indepen- 
dent element as shown in Eq. (|2l|), and its target average 
satisfies the relation 



P,,(r,r,0)(„) = (l-7V«(r))P3s(r,r,0)(„_i) + -«(r), 

(A7) 

which immediately gives the result Eq. (p2[). The same 
is obtained by setting x = x and y = y in the matrix, 



1 - 



N 



{N^ - 1) 



N 



^ l-(7V2„i)£i£) 



(A8) 



which has now the property of the probability matrix. 
The result yields 



1 ^ {N^ ~ 1){N^ - 1 + {1 - Nvir))"^) 



m \ {N^ - 1)(1 - (1 - Nv{r)y') 

(iV2 _ _ (1 _ Nv{v)Y^ 



(A9) 



APPENDIX B: MORE ON SU(2) 

As the fundamental representation of SU(2) is pseudo- 
real: t'^t°-t'^ = -vjrith Pauli matrices r", one finds 
for the Wilson lines 



'iy(x) 



.2 _ wt 



(Bl) 



Tr(W^(x)W^t(y)) ^ Xr(T2W^(x)T2T2Ty^(y)T2) 

= lTr(Ty(x)tTy(y)), (B2) 

and therefore for the eikonal factor of the singlet pair 



C4j7J(r, r, 0) = CO^(±r, , 0). 



(B3) 
(B4) 



In Ref. 1^ a model of SU(2) with the finite, fixed field 
magnitude A is studied, where the importance of the 
color degrees of freedom in the color dipole interactions 
at high energies is stressed and the spatial degrees are 
omitted. That model is hence similar to the diffraction 
case in our framework with identifying the eikonal factor 
for each zone as Wqi — cos A + 2it'^A°' sin A for the quark 
and = cos A - 2it°-A!°- sin A for the anti-quark. The 

orientation A of the field A° is assumed random and the 
average is accordingly defined by A°-J^ = \^A^^°'^ . If as- 
suming that the orientations of A and A! are uncorrelated 
and random, one finds a recursion relation 



puncor — 



= (1--Sin2 A)2puncor^^_^^ 

2 2 
-- sin^ A{\- - sin^ A). 



The solution is 



H- 2 cos 2 A 



(B5) 



(B6) 



which is the result of Ref. ||^. On the other hand, the 
case oi A = — A! where the quark and the anti-quark see 
the correlated field, which is actually simpler, one obtains 



pc 



(n) 



(1 - - sin^ 2A)Pf, 



1 



1 



1 + 2cos4A 



sin^ 2A (B7) 
(B8) 



The n- and A^-dependences of the results in both cases 
are completely consistent with Eq. (E2h with N ~ 2. 
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